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We present a systematic study of the kaon to pion multiplicity ratios {K + /tt + and K~ /rr~) in 
heavy-ion collisions from AGS to RHIC energy using the Relativistic Quantum Molecular Dynamics 
(RQMD) model. The model reasonably well describes the available experimental data on K + /tt + 
and K~ /n~ . Within the model, we find that the strong increase of the ratios with the number of 
participants is mainly due to hadronic rescattering of produced mesons with ingoing baryons and 
their resonances. The enhancement of K/tt in heavy-ion collisions with respect to elementary p + p 
interactions is larger at AGS than SPS energy, and decreases smoothly with bombarding energy. 
The total multiplicity ratios at RHIC energy are predicted by RQMD to be K + /n + = 0.19 and 
K~/ir~ = 0.15. 

PACS number(s): 25.75.-q, 24.10.Lx 
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I. INTRODUCTION 

Nuclear matter at high energy density has been ex- 
tensively studied through high energy heavy-ion colli- 
sions 0,0. The primary goal of these studies is to ob- 
serve the possible phase transition from hadronic matter 
to quark-gluon plasma (QGP) I^Q, in which quarks and 
gluons are deconfined from individual hadrons forming 
an extended region. The phase transition is predicted 
by lattice QCD calculations to occur at a temperature 
140-170 MeV and an energy density on the order of 0.5- 
1.5 GeV/fm 3 [|. It is believed that the QGP state ex- 
isted in the early universe shortly after the big bang G) , 
and may also exist in the cores of neutron stars |g| . 

If a QGP is produced in a heavy-ion collision, the 
collision system will evolve in stages from deconfined 
quarks and gluons to interacting hadrons, and finally to 
freeze-out particles which are detected. In order to ex- 
tract the information about the postulated quark-gluon 
plasma stage of heavy- ion collisions, systematic studies of 
multi-observables at freeze-out as a function of the colli- 
sion volume and bombarding energy are necessary. These 
observables include strangeness production charm 
production ||, lepton production |i"Q ], jet quenching jy], 
elliptic |Q , and other types of flow [p^ | . For a critical 
review of these observables, see Ref. U$ . 

In this article, we constrain ourselves to one of the 
above observables, namely, strangeness production. In 
particular, we focus on the kaon to pion multiplicity 
ratios (K/ir for both charge signs), since the bulk of 
strangeness produced in heavy-ion collisions are carried 
by kaons. 

It is not clear how much the K/tt ratios in heavy-ion 
reactions reflect the properties of the phase transition be- 
tween QGP and hadronic matter. Originally it was ar- 
gued that K/tt may serve as a signature of the QGP, and 



indeed early measurements at the BNL AGS and CERN 
SPS jl5| showed a significant enhancement of K + /tt + in 
heavy- ion collisions over p + p interactions. The idea was 
that, in the deconfined state, strange quark pairs (ss) 
may be copiously produced through gluon-gluon fusion 
(gg — > ss) l^yn!), while in the hadronic gas such pairs 
are produced via pairs of strange hadrons at a higher 
production threshold. Therefore the time needed for a 
hadronic gas system in kinetic equilibrium to reach chem- 
ical equilibration is significantly longer than the life time 
of a heavy-ion collision which is typically on the order 
of 10 fm [^| Jl^ , p^| . This idea, however, has been chal- 
lenged because it neglects pre-equilibrium dynamics of 
the initial stage which may considerably speed up chem- 
ical equilibration ]T^ |: Initial interactions between pro- 
duced particles and ingoing baryons are "harder" than in 
kinetic equilibrium. This is confirmed by detailed trans- 
port calculations which can reproduce the early experi- 
mental data on K+ /tt+ |l|-|l|. 

Even on the level of equilibrium physics, it is not 
entirely clear which difference to expect in a compari- 
son of the scenarios with and without phase transition. 
Strangeness in a chemically equilibrated hadronic gas 
might be as high as or higher than in a QGP |l^,|l^,^2| . 
In fact, thermal hadron gas model requires a so-called 
strangeness suppression factor to fit to experimental par- 
ticle ratios; this has been interpreted as some indica- 
tion of non-equilibrium transition from QGP to hadron 
gas [^|. In this picture, a "lower" than expected K/tt 
ratio is a signature of the QGP. 

Recently several measurements of K/tt have been made 
at both the AGS and SPS by different experiments 
p4] ps| . There has been renewed interest in the inter- 
pretation of the measured K/tt ratios ^3[|9). Interest- 
ingly, the Giessen group found the strongest deviations 
of their hadronic model (HSD) calculations from exper- 
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imental data not at SPS energy but at the lower AGS 
energy 

Despite of the theoretical mist, the K/ tt ratios may still 
be valuable observables to be studied not only to address 
questions of the phase transition but also to obtain a bet- 
ter understanding of the pre-equilibrium dynamics, the 
hadronization processes, and dynamics of hadrons in the 
medium. In these studies, the beam energy is an essential 
control parameter whose variation allows to modify two 
important variables: the initial baryon and energy den- 
sities. Both densities are of importance as to whether 
the system enters into the quark-gluon phase or remains 
hadronic all the time. If the "strangeness content" of the 
hadronic matter is substantially different from a QGP, 
discontinuity might be expected in the excitation func- 
tion of the K/tt ratios. With data from RHIC emerging 
in the near future, where QGP is likely to be formed, 
the first things to look for are changes in the observables 
departing from what have been observed at low energies. 

In this paper we present a systematic study of the 
K/tt ratios using the Relativistic Quantum Molecular 
Dynamics (RQMD) model [|30|-|32"|] . For simplicity, we 
only consider the equal and large mass nucleus-nucleus 
collisions. The goal of this paper is to provide an under- 
standing of the underlying mechanisms for K/tt enhance- 
ment by comparing model results with available experi- 
mental data at various energies. We also present predic- 
tions of the K/tt ratios at RHIC energy. Throughout the 
paper, K/tt stands for the integrated total multiplicity 
ratios unless otherwise notified. 

The paper is organized as follows. In Sec. [[J, we give 
brief description of the RQMD model. In Sec. H, we 
demonstrate that the calculated pion and kaon multi- 
plicities are in a good agreement with experimental data 
in elementary p + p interactions. In Sec. we present 
our results on K/tt ratios in calculated heavy-ion colli- 
sions. This section is divided into four parts. Part A 
presents a systematic study of K/tt as a function of the 
number of participants, parts B and C investigate the 
effects of particle rescattering and rope formation on the 
ratios, part D considers excitation functions of K/tt span- 
ning AGS, SPS, and RHIC energies (y/s « 5, 17 and 
200 AGeV, respectively), and Part E derives the K/tt 
enhancement factors. Section V summarizes our find- 
ings. 



II. THE RQMD MODEL 

RQMD is a semi-classical microscopic model which 
combines classical propagation with stochastic interac- 
tions |50|. The degree of freedom in RQMD depends 
on the relevant length and time scales of the processes 
considered. In low energy collisions (~1 AGeV), RQMD 
reduces to solving transport equations for a system of 



nucleons, other hadrons and resonances interacting in bi- 
nary collisions or via mean field. At high beam ener- 
gies (>10 AGeV), color strings and hadronic resonances 
are excited in elementary collisions; their fragmentation 
and decay lead to production of particles. Overlapping 
strings do not fragment independently but form "ropes" 
pp] , |3l| -|33|] . The secondaries which emerge from the frag- 
menting strings, ropes and resonances may re- interact. 
For all the results presented here, the so-called "cascade 
mode" is used in RQMD (no "mean field" ) . 

In order to calculate heavy-ion collisions at RHIC en- 
ergy, the RQMD code evolved from version 2.3 to 2.4. 
The physics are the same in the two versions. The change 
of version was due to technical reasons. We use version 

2.3 in calculations for all results presented in this paper, 
except for heavy-ion collisions at RHIC energy version 

2.4 is used. We have checked our AGS and SPS results 
with those from version 2.4 and our RHIC results of pe- 
ripheral heavy-ion collisions with those from version 2.3; 
consistencies were found. 

Let us note about the "rope" mechanism imple- 
mented in the model. In RQMD (20|j3|,|| (and also 
in QGSM @ and the Spanish version of DPM f§), 
strings fuse into color "ropes" if they overlap in trans- 
verse and longitudinal configuration space; the transverse 
dimension of strings is on the order of 0.8 fm ^0|. The 
strong ss enhancement in the coherent rope fields is a 
consequence of the large chromoelectric field strength, 
because quark pair production rates depend strongly on 
the ratio between the squared quark mass and the field 
strength J36|. It has been shown that the rope mecha- 
nism strongly enhances multi-strange baryon production 
p0| , p7 38 1 . However, its effect on single-strange hadron 
production is small once particle rescattering is consid- 
ered, resulting in negligible change in the K/tt ratios [j32"t . 

On the other side, rescattering between hadrons and 
resonances changes considerably the ratio between pro- 
duced strange and non-strange quarks. This change 
is about 80% for central Pb+Pb collisions at SPS en- 
ergy [p2[. Multi-step processes such as 
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are of essential importance for strangeness production, 
because intermediate resonances act as an "energy store." 
Such multi-body interactions are only frequent in a sys- 
tem of sufficient density of roughly 1 fm -3 , because 
the lifetime of the intermediate states is typically 1- 
2 fm/c pl|. Therefore, strangeness enhancement via 
rescattering is not expected to occur in the late dilute 
stages of the reactions. 
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Since feed-down affects the extracted particle yields, 
let us finally mention which hadrons are kept stable when 
we calculate the K/tt ratios. Members of the baryon 
octet (mainly A, S, Q), and the K s , K\, rj, and <f> mesons 
are not decayed in the model. In reality, all these parti- 
cles decay into pions via single and/or multi-step chan- 
nels; (f>'s and f2's decay into kaons in addition. Unless 
notified otherwise, these decay products are not included 
in the charged pion or kaon multiplicity. According to 
RQMD, the 77-decay contribution to the charged pion 
multiplicities is 2-4%; K s multiplicity is, to a good ap- 
proximation, one half of the charged kaon multiplicity, 
therefore the iCg-decay contribution to pion multiplici- 
ties can be deduced from K + /tt + and K~ /ir~; the con- 
tributions from all other sources to ir + (tt~) multiplicity 
are less than 3% (5%) and 5% (7-15%), respectively for 
p + p and heavy- ion collisions; and the 0-decay contri- 
butions are a couple to 10% in K + and 7-15% in K~ 
multiplicity in all collisions studied. 

III. RQMD RESULTS FOR P + P INTERACTIONS 

Figure [l] (top panel) shows the inclusive tt + , 7r~, K + , 
and K~ multiplicities in p + p interactions as a function 
of the center-of-mass energy (y/s) from AGS to SPS en- 
ergy. The K + /tt + and K~ /ir~ ratios are shown in the 
bottom panel. The symbols are results of RQMD calcu- 
lations. The dashed curves are parametrizations of the 
experimental p + p data f40|| , with experimental uncer- 
tainties shown in the shaded areas. Good agreement is 
found between the RQMD results and the experimental 
data, especially for yfs > 6 AGeV. We note that similar 
degrees of agreement have also been achieved for p + A 
collisions |l^^l|,^] . These agreements provide a realistic 
base for RQMD calculations of heavy-ion collisions. 

It should be noted that the agreement between the 
RQMD results and the experimental p + p data is not 
granted automatically by inputs into RQMD, such as 
parametrizations of exclusive cross sections. Depending 
on energy, resonances or strings are excited in the elemen- 
tary inelastic nucleon-nucleon collisions (see Refs. [ f32|]43"| 
for details). While the resonance parameters are taken 
from the data tables, the parameters of string fragmen- 
tation are fixed from the properties of strings created in 
e + e~ — > qq annihilations. One of the most important 
differences concerning particle production in e + e~ ver- 
sus hadron-hadron interactions is the fragmentation of 
the ingoing valence quarks (e.g., the leading-particle ef- 
fect Jll]]). It is treated in a constituent quark spectator 
fragmentation approach which keeps track of the gluonic 
"junction" connecting all three quarks of a nucleon. This 
is important for the stopping of nucleons on a nuclear 
target. 



IV. RQMD RESULTS FOR HEAVY-ION 
COLLISIONS 

We have calculated Au+Au collisions at lab beam ener- 
gies 8, 11, 14.6, 20, and 30 AGeV (-y/s = 4.3, 4.8, 5.4, 6.3, 
and 7.6 AGeV), Pb+Pb collisions at lab beam energies 
40, 60, 80, 100, 120, 140, 158, and 180 AGeV («/s = 
8.8, 10.7, 12.3, 13.8, 15.1, 16.3, 17.3, and 18.4 AGeV), and 
Au+Au collisions at RHIC energy (y/s = 200 AGeV). 
The comparisons of K + /tt + between RQMD and lower 
energy data are reported elsewhere ps)] . By studying 
these collisions, we shall try to identify the underlying 
physics within RQMD for the enhancement in K/tt ra- 
tios at the AGS and SPS energies, and present RQMD 
predictions for the ratios at RHIC energy. 

A. K/n systematics versus centrality 

In this section, we study K/tt ratios as a function of the 
collision centrality. We choose the number of participants 
(N p ) as a characterization of the collision centrality, con- 
sidering that N p is experimentally accessible via forward 
energy measurements. For RQMD results, N p is taken 
to be the number of initial nucleons that interact at least 
once with other particles. 

Preliminary experimental data on the centrality de- 
pendence of K/tt are available from AGS E866 [^4| and 
SPS NA49 ||§. Before discussing the RQMD results at 
various energies, we compare the RQMD results at AGS 
and SPS energies to the available data in Fig. |^. In order 
to facilitate the comparison, contributions from 77, <j), Kg, 
and strange (anti)baryon decays are included in the 
RQMD results. These contributions reduce the RQMD 
K + /tt + ratio by 5% (and is negligible for K~/tt~) at 
AGS energy, and by 10% for both K+/tt+ and K /tt~ 
at SPS energy. Good agreement is found between the 
RQMD results and the AGS data (within 20%). RQMD 
systematically overpredicts the ratios at SPS energy by 
10-30%. The agreement is fair with regard to the rela- 
tively large systematic errors on the data. See Parts B 
and C for further discussion. 

Figure |^ (left panel) shows K + /tt + calculated by 
RQMD as a function of N p in heavy-ion collisions at vari- 
ous bombarding energies. K + /tt + increases with y/s and 
increases from peripheral to central collisions. For all 
collision energies studied except RHIC energy, the de- 
pendence of the ratio on N p is similar. The increase in 
the value of K + /tt + is of the same magnitude from pe- 
ripheral to central collisions. Consequently, the relative 
increase in K + /tt + from peripheral to central collisions 
is larger at lower energies. 

Figure || (right panel) shows K~ /tt~ as a function of 
N p . K~/tt~ has similar dependence on N p as K + /tt + for 
collisions at energies between AGS and SPS. At each cen- 
trality, the ratio increases with </s. Unlike K + /tt + , the 
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absolute increase in K~ /ir~ from peripheral to central 
collisions is larger at higher energies. 

At RHIC energy as shown by Fig. [| the N p depen- 
dence of K + /tt + is distinctly different from those at the 
low AGS and SPS energies. Here, a saturation of the 
ratio seems to appear in central collisions. Consequently, 
the absolute increase in the ratio from peripheral to cen- 
tral collisions is smaller. On the other hand, K~ /tt~ 
continues to increase with N p , similar to those observed 
at the lower energies. The absolute value of K~ /tt~ at 
RHIC energy is significantly higher than the low energy 
values. It is also interesting to note that the absolute 
value of K + /tt + is lower at RHIC than SPS energy in 
central collisions. 

What makes K + and K~ different is, of course, the 
presence of the net baryon number. The effect of the 
fixed net baryon number is larger at lower energies. As 
a result, associate production of K + via NN — ► NK + Y 
(where Y = A or E) dominates over that of K~ via 
N N — > NK~Y. Note that pair production results in the 
same number of K +, s and K~~ 's. On the other hand, the 
difference between associate production of K + and K~ 
becomes negligible at RHIC energy; particle production 
from string and rope decays becomes important. There- 
fore, K + /n + and K~ /tt~ display a more similar shape 
and magnitude at RHIC than SPS energies, as shown in 
Fig. | 



B. Effect of particle rescattering 



It is shown in Refs. 
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within the framework of 
RQMD, that the increase in K + /tt + in heavy-ion colli- 
sions at AGS energy with respect to p + p and p + A 
is largely due to secondary rescattering among the par- 
ticles. Therefore, we shall study the K/tt ratios as a 
function of the average number of collisions (y p ) suffered 
by each nucleon. v v includes both nucleon-baryon and 
nucleon-meson interactions. Interestingly, v p shows little 
variation with bombarding energy; it depends mainly on 
the collision geometry. 

For both K + /tt + and K~/tt~, the v p dependence dis- 
plays a similar feature seen in Fig. [| These similarities 
are not surprising because the average values of v p and 
N p are correlated (although the distribution of v p ver- 
sus N p is fairly broad). For large v p , however, K + /tt + 
appears to saturate for nuclear collisions at high y/s in- 
cluding those in the SPS energy regime. 

To elaborate further the effect of particle rescattering 
on the K/ir ratios, we calculate the ratios using RQMD 
with particle rescattering turned off. To be specific, we 
include only the primary interactions of the ingoing nu- 
cleons (including rope formation) while the projectile and 
target pass through each other, but turn off the meson- 
induced interactions. The results of K + /tt + and K~ /tt~ 
are shown in the dashed curves in Fig. || for heavy-ion 



collisions at AGS and SPS energies. The RHIC energy 
results are shown in Fig. |^. For comparison, the default 
RQMD results are also shown. The top axis indicates the 
average impact parameter b corresponding to N p . The b 
values in Fig. ^ are extracted from Pb+Pb collisions at 
SPS energy, and in Fig. |^ from Au+Au collisions at RHIC 
energy. Note that the average b has a systematic devi- 
ation of ±0.5 fm for collisions at AGS energy and with 
different settings (default versus no rescattering and/or 
no rope). 

Referring to Fig. [I|and Fig. [|, the K + /tt + and K~ /tt~ 
ratios obtained with no particle rescattering show a shal- 
low increase (if not constant) with N p in heavy-ion colli- 
sions at all three energies. Thus meson-induced interac- 
tions are responsible for the significant increase in the ra- 
tios with N p obtained with the default RQMD, especially 
at AGS energy. This confirms the early results obtained 
in Refs. |^8| , ^9| . Meson- induced interactions mainly in- 
crease the kaon production rate; they do not change the 
pion multiplicity significantly. The changes in the pion 
multiplicities due to rescattering are —20%, —10-0% and 
+20% at AGS, SPS, and RHIC energies, respectively. 
Likewise, the corresponding changes in the kaon multi- 
plicity are 300%, 20-30% and 60%. Consequently, the 
K/tt ratios are increased by including meson-induced in- 
teractions. 

The above results can be understood if one keeps in 
mind that resonances play the most important role in 
order to overcome kaon production threshold. To illus- 
trate this, compare 7rA(1232) and irN collisions at the 
same relative momentum 700 MeV/c in the rest frame. 
In addition to pion production, the ttA collision has a 
certain probability to produce a KA pair. On the other 
hand, the irN collision can only produce pions. Thus ttN 
collisions do not contribute directly to the strangeness en- 
hancement. Only as a part of a many-body process do 
they play a role (e.g., by involving a g as an intermediate 
state which subsequently interacts with another hadron) . 
Similarly, the contribution from mr collisions turns out 
to be irrelevant. 

With the meson-induced interactions turned off, v p 
is significantly lower than obtained by the default 
RQMD, because only the nucleon-baryon interactions are 
counted. At AGS energy, the K/tt ratios are found to be 
constant over v p , indicating that multiple baryon-baryon 
interactions do not alter the K/tt ratios. They increase 
the individual kaon and pion multiplicities with similar 
relative magnitude. At the higher SPS and RHIC ener- 
gies, the physics in RQMD is different from that at AGS 
energy: Strings and ropes are formed and their fragmen- 
tation results in particle production. Thus, the meaning 
of v p is questionable. 

Since the rope mechanism is a novel implementation 
in RQMD in an attempt to take into account the string- 
string interactions during the hot and dense stage of 
heavy- ion collisions, we study its effect on K/tt in the 
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next section. 

C. Effect of rope formation 

The dot-dashed curves in Figs. [| and |^ show the cal- 
culated K/ir ratios in heavy-ion collisions at SPS and 
RHIC energies, respectively, by RQMD with rope for- 
mation turned off (and including particle rescattering). 
The effect of rope formation is negligible at AGS ener- 
gies. The differences between these results and the de- 
fault RQMD results (solid curves) are small. The small 
differences indicate that the effect of rope formation is 
small on K/% ratios once particle rescattering is consid- 
ered. 

It is interesting to investigate the case with both par- 
ticle rescattering and rope formation turned off, as this 
case corresponds to the simple nucleon-nucleon collisions 
with only the complication of initial multiple scattering. 
The RQMD results for this simple case are shown in the 
dotted curves in Figs. || and || for SPS and RHIC en- 
ergies, respectively. Comparing these results with the 
dashed curves (with rope and no rescattering) indicates 
that rope formation is able to make significant difference 
to K + /ir + at both SPS and RHIC energies if there is no 
subsequent rescattering. On the other hand, the effect 
on K~ /ir~ is small. 

Generally speaking, rope formation is an early phe- 
nomenon in RQMD and particle rescattering takes place 
later. Rope formation increases K + /ir + ; particle rescat- 
tering washes out these increases. Particle rescattering 
makes larger increases in K + /n + without rope formation 
(dot-dashed versus dotted curves) than with rope forma- 
tion (solid versus dashed curves). In other words, parti- 
cle rescattering and rope formation do not act on K + /n + 
additively. These observations are understandable if one 
considers rope formation a form of "string/parton rescat- 
tering" at early stages of heavy-ion collisions. 

The K~ /ir~ at SPS energy behaves differently from 
K + /tt + . Rope formation seems not to have an effect; the 
effect of particle rescattering is smaller than on K + /ir + . 
On the other hand, K~ /ir~ at RHIC energy behaves sim- 
ilarly to K + /ir + . This is expected because the net baryon 
density is low at RHIC energy. If the net baryon density 
were zero, then K + /n + and K~ /n~ would be identical. 

When both particle rescattering and rope formation 
are turned off as shown by the dotted lines in Figs. || and 
|], the K/ir ratios are rather constant over N p . Therefore, 
it is interesting to compare these results to p+p. The cal- 
culated p + p results are indicated by the shaded areas in 
Fig. | for AGS and SPS energies and in Fig. | for RHIC 
energy. Note that K + /ir + is slight lower in the isospin 
weighted nucleon-nucleon (N+N) interactions than p+p, 
while K~ /ir~ is slightly higher p5| , f44| |. As seen in the fig- 
ures, the heavy-ion results with no rescattering and no 



rope are consistent with the p+p values (with two excep- 
tions: the K+/tt+ at SPS and the at RHIC en- 
ergy) . The two exceptions are currently not understood. 
We remark that the p + p results are consistent with ex- 
trapolations of the measured heavy-ion data shown in 
Fig. | 

We now come back to the comparison between data 
and the RQMD results at SPS energy. None of the four 
curves for SPS energy shown in Fig. ^ seems to com- 
pletely agree with the data. We remark that the central 
collision data favor particle rescattering. In the following 
two sections, we only consider the default RQMD which 
include both particle rescattering and rope formation. 

D. K/n excitation functions 

We study the excitation functions of K/ir in the 
most central Au+Au/Pb+Pb collisions (N p > 350) with 
RQMD, as a necessary theoretical baseline for the search 
of abnormal behavior. The total multiplicity ratios are 
shown as open circles in Fig. || for K + /ir + (left panel) 
and K~ /w~ (right panel). The multiplicity ratios in the 
midrapidity region (1/3 < y/ybcam < 2/3) are shown as 
open squares. The midrapidity ratios are larger than the 
total multiplicity ratios. 

Data from AGS E866 and SPS NA44 @ and 

NA49 J28| are shown in the filled squares (total multi- 
plicity ratios) and filled triangles (midrapidity ratios). 
In order to exclude the effect of decay contributions in 
the data, the AGS K + /tt + data are scaled up by 5%, and 
the SPS K+/tt+ and K~ /ir~ data are scaled up by 10% 
The AGS data are well reproduced, and the SPS data are 
slightly overpredicted by RQMD. The overall agreement 
is good. 

For comparison, the results for p + p interactions from 
Fig. |l| are replotted in Fig. ||. The K/n ratios are signifi- 
cantly larger in heavy-ion collisions than the same energy 
p + p interactions. Note that the isospin weighted N + N 
interactions may be better in comparison to heavy-ion 
collisions. However, comparing heavy- ion results to p + p 
is more feasible because only p + p data are experimen- 
tally available. 

The K + /tt + ratio in central heavy-ion collisions in- 
creases at low energy, partially due to the rapid in- 
crease in kaon production near threshold in p+p interac- 
tions [|4Ci|| . However, K + /ir + saturates at high energies. 
Clearly, the continuous increase of K + /ir + with y/s in 
p + p is not seen in central heavy- ion collisions. This 
is a direct reflection of the amount of baryon stopping 
at different energies. The baryon midrapidity density 
decreases with yfs p5[ , which is well described by the 
model. The ratio at midrapidity even decreases with a/s 
at high energies. Maximum K + /tt + appears at beam 
energy ~ 40 AGeV. This energy range is currently being 
studied at the SPS. 
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As seen in Fig. ^ (right panel), K~ /ir~ continuously 
increases with y/s. The dependence is distinctly different 
from that of K + /ir + . The ratio is larger in heavy-ion 
collisions than p + p interactions at all corresponding en- 
ergies, but follow the general trend seen in p + p. The in- 
crease in K~ /tt~ from p+p to heavy-ion collisions at the 
same energy is significantly smaller than that in K + /tt + . 
The statement can be made stronger when taken into 
account that the charge asymmetry results in a higher 
K + /ir + and a lower K~/tt~ in p + p than the isospin 
weighted N + N interactions. 

At RHIC energy, respectively for K + /tt + and K~ /tt~ , 
RQMD predicts a total multiplicity ratio of 0.19 and 0.15, 
and a midrapidity multiplicity ratio of 0.19 and 0.17. The 
predictions provide a baseline for comparison to experi- 
mental data which are expected soon. 

The distinct difference between the K + /ir + and 
K~/ir~ excitation functions results from the different 
production mechanisms for K + and 1 which are con- 
nected to the presence of the baryon density |^0[ . The 
effect of the baryon density can be made more clearly in 
the kaon multiplicity ratio {K + / K~) versus y/s. Figure |^ 
shows the midrapidity ratio of K + /K~ as a function of 
•y/s. The ratio decreases steadily with ^/s. For compari- 
son, the E866 ||, NA44 @, and NA49 data ||@ are 
also shown. The data are well reproduced by RQMD. 

E. K/tt enhancement 

Now let us turn to the enhancement in K/tt ratios. We 
define the enhancement factors as the ratios of K + /ix + 
and K~ /n~ in central heavy- ion collisions over those in 
p+p interactions at the corresponding energy. Note that 
this is not a perfect definition. One may argue that a 
division by the isospin weighted N + N results is a better 
definition, but it is not experimentally preferred. 

Figure || shows the results of the enhancement in 
K + /ir + (left panel) and K~ /tt~ (right panel). The en- 
hancement is larger at low energy for both K + /n + and 
K~ /ir~ . This is partially due to the effect of kaon pro- 
duction threshold. As discussed before, meson-baryon in- 
teractions involving resonances are the primary processes 
for K/tt enhancement at AGS energy. At low energies 
near threshold, particle rescattering is most prominent. 
Take an extreme case that the beam energy is below the 
threshold: No kaons can be produced in p + p interac- 
tions; however, kaons can be produced in heavy-ion col- 
lisions by invoking the multi-step processes where more 
than two particles cooperate to overcome the production 
threshold. In this extreme case, the enhancement is in- 
finity. 

The enhancement in K + /tt + drops more quickly at 
lower t/s. Over the whole range of ^fs studied between 
AGS and SPS energies, the enhancement factor decreases 
smoothly with t/s. The quick drop at low is partially 



due to the threshold effect mentioned above. The shallow 
tail of the enhancement factor at large y/s results from 
the flattening of K + /n + at high ^fs in heavy- ion colli- 
sions and the continuous increase of the ratio with ^/s in 
p + p interactions. 

As shown in Fig. || (right panel), the enhancement in 
K~ /tt~ drops very quickly at low ^/s. The behavior is 
consistent with the high K~ production threshold. How- 
ever, at large y/s the enhancement factor is almost con- 
stant. 

Since the enhancement factors of both K + /ir + and 
K~/ir~ have a weak dependence on yfs at high ener- 
gies, it is possible, without RQMD calculations for RHIC 
energy, to predict the enhancement factors at RHIC en- 
ergy by extrapolation. In order to do so, we parametrize 
the enhancement factor as a function of y/s. We choose 
to parametrize the K + /tt + enhancement factor by the 
functional form A + B/y/s. Fits to the enhancement 
factors shown in the left panel of Fig. || (of course ex- 
cluding that at RHIC energy) yield 1.8 + 11.2GeV/^/s 
and 1.8 + 16.8GeV/y / s, respectively for the enhancement 
in the overall K + /ir + and the midrapidity K + /tt + . The 
X 2 /Ndf of the fits are on the order of 1.5. The fit re- 
sults are superimposed in Fig. || (left panel) as the dashed 
and dotted curves. The two curves essentially overlap at 
RHIC energy, and slightly miss the calculated results. 

Because of the sharp drop of the K~/tt~ enhance- 
ment factor at low energy, we choose to parametrize 
the K~ /t:~ enhancement factor by the functional form 
A + S/(Vs - 2.87GeV) 2 , where 2.87 GeV is the kaon 
pair production threshold in NN — > NNK + K~ . Fits to 
the enhancement factor in the right panel of Fig. ^| (ex- 
cluding that at RHIC energy) yield 1.4 + 6.4GeV/(Vs - 
2.87GeV) 2 and 1.7 + 9.1GeV/(Vs - 2.87GeV) 2 , respec- 
tively for the enhancement in the overall K~ /tt~~ and the 
midrapidity K~ /tt~. The x 2 /Ndf of the fits are on the 
order of 2.5. The large x 2 /Ndf is mainly due to the 
lowest y/s point. If the lowest y/s point is excluded from 
the fits, x 2 /Ndf becomes about 1.5 without essential 
change in the fit results. The fit results are superim- 
posed in Fig. H (right panel) as the dashed and dotted 
curves. They can satisfactorily describe the enhancement 
factor at RHIC energy as directly calculated by RQMD. 

V. CONCLUSIONS 

We have reported a systematic study of the kaon to 
pion multiplicity ratios, K + /ti + and K~ /tt~, in heavy- 
ion collisions as a function of the bombarding energy from 
AGS to RHIC energy, using the RQMD model. We have 
demonstrated that the kaon and pion multiplicity data in 
p+p interactions are well reproduced by RQMD, validat- 
ing the comparison between data and RQMD for heavy- 
ion collisions. 
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The RQMD-calculated ratios in heavy-ion collisions 
are higher than those in p + p interactions at the same 
energy, and increase from peripheral to central collisions. 
By comparing the results to those calculated by RQMD 
with particle rescattering (meson-induced interactions) 
turned off, we conclude that the K/tt enhancement in 
central collisions with respect to peripheral collisions and 
p + p interactions is mainly due to meson-induced inter- 
actions, especially at the low AGS energy. We further 
conclude that rope formation does not change the K/tt 
ratios significantly once particle rescattering is consid- 
ered. 

It is found that the K/tt enhancement in central heavy- 
ion collisions over p+p interactions is larger at AGS than 
SPS energy, and decreases smoothly with bombarding en- 
ergy. This behavior is consistent with the combination 
of the threshold effect of kaon production in p + p in- 
teractions and the dropping baryon density in heavy-ion 
collisions with increasing bombarding energy. 

The RQMD model reasonably describes the available 
K + /ir + and K~ /tt~ data in heavy-ion collisions. The 
observed K/tt enhancement at the AGS and SPS can be 
understood in the RQMD model with hadronic rescat- 
tering and string degrees of freedom. 

The RQMD K/tt results agree with experimental data 
better at midrapidity in the most central collisions. This 
is not surprising since it is the case where equilibrium 
is most likely to be reached. There the details of the 
description in any model is less relevant. 

The total multiplicity ratios at RHIC energy are pre- 
dicted by RQMD to be K+ /tt+ = 0.19 and K~ /tt~ = 
0.15. The midrapidity ratios at RHIC energy are pre- 
dicted by RQMD to be K+ /tt+ = 0.19 and K~ /tt~ = 
0.17. 
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FIG. 1. RQMD results of p+p interactions at various cen- 
ter-of-mass energies (y/s). Top panel: ir + (filled circles), ir~ 
(open circles), K + (filled squares), and K~ (open squares) 
total multiplicities. Bottom panel: K + (filled circles) 
and K~ /n~ (open circles) ratios. All statistical errors on 
the RQMD results are smaller than the symbol size. The 
dashed curves are parametrizations of experimental data [^oj . 
The shaded areas around the curves in the top panel reflect 
15% experimental uncertainties. Those in the bottom panel 
reflect the propagated 21% uncertainties (only one half of the 
experimental uncertainties on K + /tt + is shown for clarity). 
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FIG. 2. Comparison of the kaon to pion total multiplicity ratios, K + /ir + (left panel) and K~ /n~ (right panel), between 
RQMD results and preliminary data from Au+Au/Pb+Pb collisions at AGS and SPS energies as a function of the number of 
participants (N p ). The RQMD results are shown in the small symbols connected by the lines. The filled circles are preliminary 
data from AGS E866 p4| . The open squares are preliminary data from SPS NA49 p^ |. In order to facilitate the comparison, 
the RQMD results are scaled down to include decay contributions. Errors (either shown or smaller than the symbol size) are 
statistical only for both the RQMD results and the E866 data. Errors on the NA49 data are quadratic sum of statistical and 
5% systematic errors. 
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FIG. 3. RQMD results of the kaon to pion total multiplicity ratios, K + /ir + (left panel) and K~ /tt~ (right panel), in 
Au+Au/Pb+Pb collisions at various bombarding energies as a function of the number of participants (JV P ). The lines are 
labeled by the lab beam energies except for RHIC where y/s = 200 AGeV. Errors are statistical only, either shown or smaller 
than the symbol size. The top axis indicates the average impact parameter (b) for Au+Au collisions at RHIC energy from 
default RQMD. The average impact parameters for other energies have a systematic variation of ±0.5 fm from these values. 
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FIG. 4. RQMD results of the kaon to pion total multiplicity ratios, K + /-k + (left panel) and K~ /ir~ (right panel), in Au+Au 
collisions at AGS energy and in Pb+Pb collisions at SPS energy as a function of the number of participants (N p ). Shown 
are the RQMD results with meson-induced interactions and rope formation both turned off (dotted line), with meson-induced 
interactions turned off but rope formation included (dashed lines), with meson- induced interactions included but rope formation 
turned off (dot-dashed line), and with the default setting (solid lines). The AGS results are shown with the open circles, and 
the SPS results are shown with the solid squares. Rope formation has negligible effect at AGS energy. Errors are statistical 
only, either shown or smaller than the symbol size. The shaded areas indicate the corresponding ratios in p + p interactions 
calculated by default RQMD at AGS (lower) and SPS (upper) energies. The top axis indicates the average impact parameter 
(6) for Pb+Pb collisions at SPS energy from default RQMD. The average impact parameters for AGS energy and for the other 
settings have a systematic variation of ±0.5 fm from these values. 
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FIG. 5. RQMD results of the kaon to pion total multiplicity ratios, K + /tt + (left panel) and K~ (right panel), in Au+Au 
collisions at RHIC energy as a function of the number of participants (N p ). Shown are the RQMD results with meson-induced 
interactions and rope formation both turned off (dotted line), with meson- induced interactions turned off but rope formation 
included (dashed line), with meson-induced interactions included but rope formation turned off (dot-dashed line), and with 
the default setting (solid line). Errors are statistical only, either shown or smaller than the symbol size. The shaded areas 
indicate the corresponding ratios in p + p interactions calculated by default RQMD. The top axis indicates the average impact 
parameter (b) for Au+Au collisions at RHIC energy from default RQMD. The average impact parameters for the other settings 
have a systematic variation of ±0.5 fm from these values. 



0.3 



0.25 



0.2 



+ ^ 0.15 



0.1 



0.05 



AGS 



SPS 

I 



RHIC 



A AA midrapidity (data) 

■ AA overall (data) 

□ AA midrapidity (RQMD) - 

O AA overall (RQMD) 

p+p (RQMD) 



0.2 



10 



10 



[GeV] 




0.02 : oft 



10 



RHIC 

I 

□ 

o 



A AA midrapidity (data) 

■ AA overall (data) 

□ AA midrapidity (RQMD) 

O AA overall (RQMD) 

p+p (RQMD) 



10 



Vs [GeV] 



FIG. 6. RQMD results of the kaon to pion multiplicity ratios, K + /-k + (left panel) and K~ /n~ (right panel), in 
Au+Au/Pb+Pb central collisions as a function of the nucleon-nucleon center-of-mass energy (y/s). The total multiplicity 
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data {mJ. Data from AGS E866 J||§|] and SPS NA44 §| and NA49 §| are also shown for the total (filled square) and 
midrapidity (filled triangle) multiplicity ratios. These data have been scaled up to approximately exclude decay contributions. 
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than the symbol size) are statistical only for both the RQMD results and the E866 data. Errors on the NA44 data are combined 
statistical and systematic errors. Errors on the NA49 data are quadratic sum of statistical and 5% systematic errors. 
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FIG. 8. RQMD results of the strangeness enhancement factor, the ratios of K + /n + (left panel) and K~ /ir~ (right panel) in 
central heavy-ion collisions over those for the same energy p+p interactions, as a function of the nucleon-nucleon center-of-mass 
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